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Foreword 


The  Navy’s  Fleet  Numerical  Oceanography  Center  (FNOC)  in  Monterey, 
California,  is  responsible  for  providing  three-dimensional  global  and  regional 
scale  analyses  of  ocean  temperature  to  the  Fleet.  These  products  have  a  direct 
and  important  impact  on  the  Navy’s  mission  and  daily  operations.  However, 
the  amount  of  oceanographic  data  that  exists  is  limited,  at  best,  due  to  the 
vast  ocean  area  that  must  be  sampled  and  the  relatively  small  number  of  ships 
and  buoys  available  to  make  observations.  While  remotely  sensed  measure¬ 
ments  are  more  abundant,  large  areas  of  the  oceans  are  normally  obscured 
by  cloud  cover.  Furthermore,  the  only  satellite  data  being  directly  used  for 
temperature  analysis  at  this  time  are  the  Multichannel  Sea  Surface 
Temperatures,  which  provide  information  only  about  the  ocean’s  surface. 
Because  of  the  data  sparsity,  especially  below  the  surface,  it  is  essential  that 
we  make  optimum  use  of  each  observation. 

The  Naval  Ocean  Research  and  Development  Activity  has  worked  closely 
with  scientists  at  FNOC  to  develop  a  new  and  improved  thermal  analysis 
product,  the  Optimum  Thermal  Interpolation  System  (OTIS),  to  meet  this 
goal.  This  report  describes  in  detail  the  state-of-the-art  techniques  that  are 
employed  in  OTIS  to  produce  a  nowcast  of  the  three-dimensional  global  ocean 
thermal  structure.  While  OTIS  is  a  useful  environmental  product  in  its  own 
right,  the  output  from  OTIS  provides  the  FNOC  acoustic  models  with  improved 
input  data  fields.  Thus,  the  implementation  of  OTIS  as  FNOC’s  operational 
analysis  has  resulted  in  the  delivery  of  superior  environmental  and  acoustic 
products  to  the  Fleet. 


W.  B.  Moseley 
Technical  Director 


J.  B.  Tupaz,  Captain,  USN 
Commanding  Officer 


Executive  Summary 


The  Optimum  Thermal  Interpolation  System  (OTIS)  is  an  ocean  thermal 
analysis  product  developed  for  real-time  operational  use  at  the  Fleet  Numerical 
Oceanography  Center.  OTIS  is  expected  to  become  the  centerpiece  of  the 
Navy’s  ocean  thermal  analysis  and  prediction  capabilities  both  ashore  and 
afloat.  It  provides  a  rigorous  framework  for  the  synergistic  combination  of 
real-time  data,  climatology,  and  predictions  from  ocean  mixed-layer  and 
circulation  models  to  produce  the  Navy’s  most  accurate  representation  of  ocean 
thermal  structure  on  global  and  regional  scales.  OTIS  is  particularly  well  suited 
for  utilization  of  remotely  sensed  data  from  satellites  because  of  its  ability 
to  account  for  the  relative  accuracies  of  various  types  of  data. 

OTIS  is  based  on  the  optimum  interpolation  (OI)  data  assimilation 
methodology.  Basically,  the  OI  technique  maps  observations  distributed 
nonuniformly  in  space  and  time  to  a  uniformly  gridded  synoptic  representation, 
or  analysis,  of  the  target  field.  The  analysis  is  constructed  as  a  first-guess 
background  field  plus  an  anomaly  relative  to  that  field.  The  analyzed  anomaly 
at  a  particular  gridpoint  is  given  by  a  weighted  combination  of  observed  and 
model-predicted  anomalies,  with  the  space-time  autocorrelation  function  for 
the  resolvable  anomalies  governing  which  observations  contribute.  The  OI 
technique  provides  the  optimum  weights  applied  to  each  anomaly,  such  that 
the  resulting  analysis  error  will  be  minimized  in  a  least-squares  sense.  The 
technique  also  provides  an  estimate  of  this  error. 

The  basic  inputs  to  OTIS  are  the  statistics  defining  both  the  resolvable  and 
subgrid-scale  variability  of  the  target  field  about  the  background  field, 
the  instrumental  error  characteristics  of  the  measurement  systems  providing  the 
observations,  and  the  error  characteristics  of  the  forecast  model  that  supplies 
the  predictions.  The  end  product  is  a  statistically  optimum  gridded  represen¬ 
tation  of  the  current  ocean  thermal  structure,  which  may  in  turn  be  used 
to  initialize  prognostic  thermodynamic  and  ocean  circulation  models,  or  to 
compute  sound-speed  profiles  for  input  to  ocean  acoustic  models. 
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Surface  infrared  radiative  heat  flux. 

Atmospheric  drag  coefficient. 

Specific  heat  of  seawater. 

Damping  coefficient  for  inertial  oscillations. 

Downward  flux  of  solar  radiation. 

Solar  radiation  penetrating  the  sea  surface. 

Coriolis  parameter. 

Surface  sensible  heat  flux. 

Surface  latent  heat  flux. 

Salinity. 

Salinity  at  the  sea  surface. 

Temperature. 

Reference  climatological  SST  gradient. 

Analyzed  temperature  at  gridpoint  k. 

Climatological  temperature  at  gridpoint  k. 

Observed  temperature  at  location  /. 

TOPS-predicted  temperature  at  gridpoint  k. 

Temperature  resulting  from  previous  OTIS  analysis. 

Individual  MCSST  observation  at  location  j. 

Weight  given  to  observation  /  assimilated  at  gridpoint  k. 

Weight  given  to  TOPS  prediction  at  gridpoint  k. 

Weight  given  to  climatology  at  gridpoint  k. 

Background  eddy  diffusion  coefficient. 

Number  of  observations  assimilated  at  gridpoint  k. 

Correlation  between  resolvable  thermal  anomalies  at  location 
and  location  j  in  space-time. 

Noise-to-signal  ratio  for  observation  /. 

Noise-to-signal  ratio  for  TOPS  prediction  at  gridpoint  k. 

Noise-to-signal  ratio  for  subgrid-scale  error  at  location  i. 

Average  of  A"  over  entire  horizontal  analysis  domain. 

Horizontal  eddy  diffusion  coefficient. 


p  Reference  density  for  seawater. 

pa  Reference  density  for  air  near  the  surface. 

Ak  Inverse  of  east-west  decorrelation  scale  of  resolvable  thermal 

anomalies  at  gridpoint  k. 

Bk  Inverse  of  north-south  decorrelation  scale  of  resolvable 

thermal  anomalies  at  gridpoint  k. 

Ck  Inverse  of  temporal  decorrelation  scale  of  resolvable  thermal 

anomalies  at  gridpoint  k. 

D  Number  of  analysis  cycles  (days)  for  temperature  anomaly 

to  decay  to  10°7o  of  initial  value  at  a  data  void  gridpoint. 

U,  V  x-  and  y-components  of  wind  vector  at  height  of  19.5  m 

above  the  sea  surface. 

u,  v  x-  and  y-components  of  drift  current. 

ug,vg  x-  and  y-components  of  the  geostrophic  current. 

ua,va  x-  and  y-components  of  the  advection  current. 

wa  z-component  of  the  advection  current. 

LX  Reference  grid  space. 

Ax,y  East-west  distance  between  location  /'  and  location  j. 

A ytJ  North-south  distance  between  location  /  and  location  j. 

A t'j  Time  difference  between  observation  i  and  observation  j. 

(A7-*  )„  Change  in  TOPS  predicted  temperature  from  level  n- 1  to 
level  n  (level  0  implies  the  surface)  at  gridpoint  k. 

A TClm  Climatological  change  in  temperature  between  previous  and 
current  analysis  times. 

hk  Analyzed  mixed-layer  depth  at  gridpoint  k. 

hk  Climatological  mixed-layer  depth  at  gridpoint  k. 

h°  Observed  mixed-layer  depth  at  gridpoint  k. 

o°  RMS  observational  error  of  observation  /. 

ok  RMS  error  of  TOPS  prediction  at  gridpoint  k. 

o‘  Long-term  RMS  variation  of  resolvable  thermal  anomalies 

about  climatology  at  location  /. 

o,ft  RMS  difference  between  bathy  observations  and  climatology 

at  location  i. 


v 


o  '  RMS  difference  between  bathy  observations  and  climatology 

calculated  over  entire  horizontal  analysis  domain. 

o'k  Expected  RMS  temperature  difference  between  Tk  and 

independent  bathy  observations. 

o,"  RMS  subgrid  scale  error  at  location  /. 

o'  RMS  error  added  to  MCSST  super  ob  as  a  result  of  the  way  in 

which  the  block  averaging  of  the  MCSST  data  is  done. 

o[  RMS  error  of  resolvable  thermal  field  in  resulting  analysis  at 

gridpoint  k. 

o'k  Steady-state  value  for  o°k  at  gridpoint  which  is  devoid  of  data  for 
a  long  period  of  time. 

( o  )tx  RMS  error  added  to  /th  bathy  observation  below  depth  zB 
by  the  downward  extrapolation  of  the  bathy. 

v.  Anomaly  reduction  factor  at  gridpoint  k. 

M  Number  of  MCSST  observations  averaged  around  a  gridpoint 

to  form  a  “super  ob.  ' 

£  Standard  deviation  of  instrumental  error  for  MCSST,  SHIP,  and 

BATHY  observations. 

c  Mean  of  instrumental  error  (i.e.,  bias)  for  MCSST  observations. 

(/,  RMS  TOPS  error  growth  over  one  forecast  cycle  at  gridpoint  k. 

Depth,  measured  positive  downward  from  the  surface  in  meters. 
Deepest  reported  depth  in  a  bathy  observation. 
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Technical  Description  of  the  Optimum  Thermal  Interpolation  System, 
Version  1:  A  Model  for  Oceanographic  Data  Assimilation 


1.0  Introduction 

The  Optimum  Thermal  Interpolation  System  (OTIS) 
is  an  oeean  thermal  analysis  product  developed  for  real¬ 
time  operational  use  at  the  Fleet  Numerical 
Oceanography  Center  (FNOC).  The  development  of 
OTIS  Version  1  is  described  by  Innis  (1983),  Innis 
(1985),  Innis  and  Mendenhall  (1986),  and  STX  (1987); 
the  OTIS  1  software  is  described  by  SAIC  (1988a), 
SAIC  (1988b),  and  SAIC  (1988c). 

OTIS  is  expected  to  become  the  centerpiece  of  the 
Navy’s  ocean  thermal  analysis  and  prediction 
capabilities  both  ashore  and  afloat.  It  provides  a 
rigorous  framework  for  the  synergistic  combination 
of  real-time  data,  climatology,  and  predictions  from 
ocean  mixed-layer  and  circulation  models  to  produce 
the  Navy’s  most  accurate  representation  of  ocean 
thermal  structure  on  global  and  regional  scales.  OTIS 
is  particularly  well  suited  for  utilization  of  remotely 
sensed  data  from  satellites.  It  can  run  on  coarse- 
resolution  grids  with  global  coverage  or  on 
eddy-resolving  grids  with  regional  coverage.  The 
OTIS  1  software,  representing  eight  main  programs 
and  approximately  20,000  lines  of  source  code,  is  well- 
documented  and  fully  integrated  into  the  FNOC 
computer  system,  and  is  under  configuration 
management  control.  This  report  provides  a  complete 
and  detailed  technical  description  of  OTIS  1  as  it  is 
now  implemented  at  FNOC. 

OTIS  is  based  on  the  Optimum  Interpolation  (Ol) 
data  assimilation  technique  of  Gandin  (1965),  which 
is  used  widely  in  meteorology.  The  OI  technique  has 
been  applied  to  oceanographic  analysis  by  Bretherton 
et  al.  (1976),  Freeland  and  Gould  (1976),  McWilliams 
(1976),  White  (1977),  Dorman  and  Saur  (1978),  White 
and  Bernstein  (1979),  Bernstein  and  White  (1981), 
Carter  and  Robinson  (1981),  Clancy  (1983),  Robinson 
(1983),  Roemmich  (1983),  Robinson  and  Leslie  (1985), 
McWilliams  et  al.  (1986),  Hua  et  al.  (1986),  Robinson 
(1986),  Robinson  and  Walstad  (1986),  Robinson  et  al. 
(1986),  Hansen  and  Herman  (1986),  Cummings  (1986), 
Carter  and  Robinson  (1987),  Beasley  (1987),  Carton 
(1987),  Murray  (1987),  Robinson  et  al.  (1987),  Carter 
et  al.  (1988),  Grillak i-Steiert  and  Amoroso  (1988),  and 
others. 


Basically,  the  OI  technique  maps  observations 
distributed  nonuniformly  in  space  and  time  to  a 
uniformly  gridded  synoptic  representation,  or 
analysis,  of  the  target  field.  As  in  most  objective  analy¬ 
sis  techniques,  the  concept  of  resolvable  versus 
subgrid-scale  features  is  fundamental  in  OI.  Resolvable 
features  are  those  of  spatial  extent  greater  than  twice 
the  mesh  length  of  the  analysis  grid,  and  thus  are 
capable  of  being  at  least  marginally  represented  in 
discrete  fashion  by  the  grid.  Subgrid-scale  features  are 
those  with  spatial  extent  less  than  twice  the  mesh  length, 
which,  as  a  result,  cannot  be  resolved  by  the  grid.  As 
far  as  the  analysis  is  concerned,  subgrid-scale  features 
contribute  to  the  error  in  the  observations. 

An  OI  analysis  is  constructed  as  a  background  or 
“first-guess”  Field  plus  an  anomaly  relative  to  this  field. 
The  anomaly  at  a  particular  gridpoint  is  given  by  a 
weighted  combination  of  observed  and  model-predicted 
anomalies,  with  the  space-time  autocorrelation 
function  for  the  resolvable  anomalies  governing  which 
observations  contribute.  The  OI  technique  provides  the 
optimum  weights  applied  to  the  observations  and 
predictions,  such  that  the  resulting  analysis  error  will 
be  minimized  in  a  least-squares  sense.  The  technique 
also  provides  an  estimate  of  this  error.  The  basic  inputs 
to  this  process  are  the  statistics  defining  both  the 
resolvable  and  subgrid-scale  variability  of  the  target 
field  about  the  background  field,  the  instrumental  error 
characteristics  of  the  measurement  systems  providing 
the  observations,  and  the  error  characteristics  of  the 
forecast  model  that  supplies  the  predictions. 

The  remainder  of  this  report  is  organized  as  follows. 
Section  2  gives  a  technically  detailed  description  of  the 
data  assimilation  technique  used  in  OTIS.  Section  3 
discusses  the  user-specified  free  parameters  required 
by  the  technique.  Section  4  presents  the  data-screening 
algorithm  used  by  OTIS  to  filter  out  spurious  data. 
Section  5  provides  a  technical  overview  of  the  ocean 
prediction  models  used  in  conjunction  with  OTIS. 
Section  6  presents  a  series  of  idealized  experiments  that 
provide  insight  into  how  OTIS  chooses  the  relative 
weighting  of  climatology,  data,  and  model  prediction 
to  arrive  at  the  analyzed  temperature.  Section  7 
discusses  the  horizontal  grids  on  which  OTIS  can 
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function  and  Section  8  discusses  the  climatologies  it  can 
use.  Section  9  is  a  brief  overview  of  the  use  of  “ocean 
feature  models,”  which  allow  OTIS  to  map  the  subsur¬ 
face  structure  of  fronts  and  eddies  accurately  from  only 
their  surface  signatures.  Section  10  shows  examples  of 
output  generated  by  OTIS.  Section  II  discusses  the 
ongoing  and  automatic  validation  of  OTIS  against 
independent  bathythermograph  (hereafter  referred  to 
as  “bathy”)  data  at  FNOC.  Finally,  Section  12  summa¬ 
rizes  the  central  role  OTIS  plays  in  linking  other  data 
bases  and  models  and  discusses  future  developments. 

Readers  unconcerned  with  the  technical  details  of 
OTIS  should  skip  to  Section  7. 

2.0  Description  of  the  Analysis  Procedure 

OTIS  can  be  executed  in  either  “shallow”  or  “deep” 
analysis  modes.  In  the  shallow  mode,  the  deepest 
analyzed  depth  is  400  m,  while  in  the  deep  mode  it  is 
5000  m.  In  performing  an  analysis,  OTIS  proceeds 
point-by-point  through  its  horizontal  grid,  producing 
a  complete  profile  from  the  surface  to  the  deepest 
analyzed  depth  at  each  gridpoint  before  advancing  to 
the  next.  At  each  of  these  gridpoints,  the  analysis 
sequent  is  (1)  sea  surface  temperature  (SST), 
(2)  temperature  in  the  mixed  layer,  and  (3)  temperature 
below  the  mixed  layer. 


optimum  analysis  (i.e.,  one  that  will  have  a  minimal 
mean-square  error  when  averaged  over  a  large  number 
of  analysis  realizations)  are  obtained  by  solving  the 
following  set  of  N+  1  equations 

N 

*i*,  =  I  (n ,j  +  6,jK')akJ  +  (n,*  +  6ikK)Pk 

j  =  i 

(2-2) 

for  i  =  1,2, .. .  N,(N  +  1)  =  A , 

where  rjy  is  the  space-time  autocorrelation  between 
climatological  anomalies  at  location  /  and  location  j 
(“location  A ”  implies  the  gridpoint).  A"  is  the 
noise-to-signal  ratio  for  the  /th  observation,  Xk  is 
the  noise-to-signal  ratio  for  the  TOPS  prediction,  and 
6jj  is  the  Kronecker  delta  function  defined  as 

6,,  =  1  for  /  =  j 
6,j  =  0  for  /  #  y . 

The  autocorrelation  function  r\tJ  is  taken  to  be 
r) ,j  =  exp {-A\  A xfj  -  B\ by*  -  C\ A / j )  (2-3) 


2.1  Sea  Surface  Temperature  Analysis 

OTIS  follows  the  approach  of  Alaka  and  Elvander 
(1972)  by  using  climatology  as  the  first-guess  field  and 
treating  input  from  a  model  prediction  as  a  special  type 
of  observation.  The  model  prediction  is  normally 
produced  by  the  Thermodynamic  Ocean  Prediction 
System  (TOPS)  mixed-layer  model  described  in 
Section  5.  Thus,  OTIS  represents  the  analyzed  SST 
at  the  Ath  gridpoint  Tk  as 

r:  -  rk  +  £  akl(7?  -  r‘)  +  pk(Tkp  -  rk) , 

(2-1) 

where  J'l  is  the  climatological  temperature  at  the 
gridpoint  valid  at  analysis  time,  Tk  the  predicted 
temperature  at  the  gridpoint  produced  by  a  24-hour 
TOPS  forecast  from  the  previous  day’s  OTIS  analy¬ 
sis,  77'the  observed  temperature  at  location  i,  T- 
the  eiimatological  temperature  at  location  /  valid  at  the 
time  of  observation  /,  akj  the  weight  applied  to  the  /th 
observed  climatological  anomaly  assimilated  at  the 
gridpoint,  fik  the  weight  applied  to  the  TOPS- 
prcdicted  anomaly  at  the  gridpoint,  and  N  the  number 
of  observations  assimilated  at  the  gridpoint. 

From  the  basic  theory  of  OI  (e.g.,  see  Alaka  and 
Elvander,  1972,  or  the  Appendix  of  this  report),  the 
weights  akl  and  fik  necessary  to  produce  a  statistically 


where  Ax,y,  A ytj  and  A ttJ  are  the  east-west,  north- 
south,  and  time  separations  between  locations  /  and 
j  in  space-time;  Ak,  Bk  and  Ck  are  the  inverses  of  the 
east-west,  north-south  and  time  decorrelation  scales; 
and  location  k  implies  the  gridpoint  at  analysis  time. 
This  functional  form  for  rj  has  been  used  widely  in 
meteorological  and  oceanographic  analysis.  The 
quantities  A,  B  and  C  are  user-specified  parameters 
that  must  be  chosen  to  reflect  the  ocean’s  resolvable 
variability  about  climatology.  They  tend  to  control  the 
space-time  domain  over  which  an  individual 
observation  will  influence  the  analysis. 

The  quantities  A"  and  A£  in  equation  (2-2)  are 
given  by 


and 


K 


(2-5) 


where  0°  is  the  root-mean-square  (RMS)  error  of  the 
/th  observation,  ak  is  the  RMS  error  of  the  TOPS 
prediction  at  gridpoint  A,  and  o,fand  ock  are  the  RMS 
variations  of  the  thermal  field  about  climatology  at 
location  /  and  gridpoint  A,  respectively. 


Determination  of  o"  for  use  in  equation  (2-4) 
depends  primarily  on  what  type  of  data  the  /th  obser¬ 
vation  is.  Three  types  of  data  contribute  to  the  OTIS 
SST  analysis:  ship  injection/bucket  temperature 
measurements,  bathythermograph  measurements,  and 
Multi-Channel  Sea  Surface  Temperature  (MCSST) 
measurements  from  the  Advanced  Very  High  Resolu¬ 
tion  Radiometer  (AVHRR)  instrument  aboard  the 
National  Oceanic  and  Atmospheric  Administration’s 
polar-orbiting  satellites.  For  ship  and  bathy  measure¬ 
ments,  o''  is  given  by 


(°;  htiip  £.s ////■  +  (°<  )  (2-6) 

and 

(°/  hiuin  ~  f-/;  > my  +  (°i  )  (2-7) 


where  f  Sl!ll,  is  the  standard  deviation  of  the 
instrumental  error  foi  the  ship  SST  measurements, 
f-r ,.Ull)  is  the  standard  deviation  of  the  instrumental 
error  for  the  bathy  measurements,  and  o"  is  the  RMS 
subgrid-scale  error  at  location  /  (i.e.,  the  RMS  value 
of  thermal  anomalies  with  spatial  extent  less  than  twice 
the  mesh  length  of  the  analysis  grid  and  thus 
unresolvable  by  the  grid).  The  quantities  tSHIP  and 
iBArHY  are  user-specified  parameters  that  include  the 
effect  of  basic  instrumental  error,  as  well  as  report  and 
transmit  errors.  Note  that  the  ship  and  bathy  errors 
are  assumed  to  have  negligible  means,  so  that  iSH,P 
and  zBA  fH)  can  also  be  thought  of  as  the  RMS  errors. 
The  values  of  cSflIP  and  £BATHY  are  chosen  to  reflect 
the  general  reliability  of  these  data  and  tend 
to  control  their  weighting  relative  to  climatology,  other 
data,  and  TOPS  in  the  analysis. 

The  MCSST  data  are  handled  in  a  special  way,  which 
complicates  the  determination  of  (o°)^CSS7.  Before 
assimilation  into  the  analysis,  the  MCSST  reports 
within  one-half  mesh  length  of  each  gridpoint  are  block 
averaged  around  the  gridpoint  to  form  “super  obs” 
assumed  valid  at  the  gridpoint.  This  averaging  process 
is  done  because  the  high  MCSST  data  densities  usually 
result  in  many  MCSST  reports  around  a  gridpoint,  and 
this  process  provides  an  efficient  way  to  incorporate 
these  data  without  losing  any  information  on  the 
resolvable  thermal  field.  Thus,  if  the  /th  observation 
in  equation  (2-1)  is  an  MCSST  observation,  then  it  is 
in  fact  given  by 


l  v' 

( T,  ft/r.vsr  ~  m  I  T;  - 


(2-8) 


where  t 1  is  the  yth  individual  MCSST  report  and  M 
is  the  total  number  of  individual  MCSST  reports  used 
in  the  block  average.  Thus,  in  view  of  equation  (2-8), 
the  expression  for  (o")2uc SST  becomes 


(°;  ) MCSST  ~~  foj  [£  MCSST  +  (°i  )  +  (°,  )  1  +  e.\ 


\nssi 


(2-9) 


where  tMCSSI  and  e%tcssr  are  the  standard  deviation 
and  mean,  respectively,  of  the  instrumental  error  for 
the  MCSST  reports,  and  o'  is  the  standard  deviation 
of  the  SST  field  resulting  from  resolvable  SST 
variations  over  the  domain  in  which  the  MCSST  block 
averaging  is  done.  The  term  (o')2  represents  the 
variance  added  because  the  differencing  of  the  MCSST 
data  with  climatology  is  done  after  the  block  averaging 
and  with  the  climatological  temperature  valid  at  the 
gridpoint,  rather  than  before  the  block  averaging  and 
with  temperatures  obtained  by  interpolating 
climatology  to  each  individual  MCSST  observation. 
This  term  is  estimated  by  assuming  a  constant  reference 
SST  gradient  and  averaging  the  resulting  temperature 
variance  over  one  grid  space.  Thus, 


(°,T 


(T*AX)2 
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(2-10) 


where  7"*  is  a  reference  climatological  SST  gradient 
(taken  to  be  0.01  °C  km'1 )  and  AX  is  a  reference  grid 
space  for  the  OTIS  grid.  Thus,  for  the  FNOC  63  x  63 
hemispheric  grid  (AX  =  320  km),  o'  =  0.85°C. 

No  attempt  is  made  to  correct  for  any  mean  error 
(i.e.,  bias)  in  the  MCSST  data.  Rather,  eAICSST 
effectively  becomes  a  lower  bound  on  the  MCSST 
observational  error  when  M  becomes  large.  The 
quantities  £A/cssr  and  eMCSST  are  user-specified 
parameters  that  tend  to  control  the  relative  weighting 
of  the  MCSST  data  in  the  analysis. 

The  individual  MCSST  reports  that  contribute  to  the 
super  obs  have  varying  ages,  extending  up  to  the  length 
of  time  defining  the  temporal  window  for  this  type  of 
data  (see  Section  4).  The  age  assigned  to  the  resulting 
super  obs,  however,  is  set  to  48  hours  relative  to  the 
analysis  date-time  group. 

The  value  of  ok  for  use  in  equation  (2-5)  is 
obtained  from 


( °k )  ~  ( °k\rvmAL  +  Gk  .  (2-11) 

where  (ok),SITIAL  is  the  RMS  error  of  the  resolvable 
thermal  field  from  the  previous  analysis  initializing 
TOPS,  and  Gk  is  the  RMS  TOPS  error  growth  over 


one  analysis-forecast-analysis  cycle  (i.e.,  24  hours). 
Thus,  following  the  approach  of  Bengtsson  and 
Gustafsson  (1972),  the  error  in  the  predicted  thermal 
field  is  taken  to  be  the  error  in  the  initial  conditions 
plus  the  error  added  due  to  error  growth  in  the  forecast 
model.  The  quantity  G  is  a  user-specified  parameter 
that  must  be  chosen  to  reflect  the  rate  of  error  growth 
in  TOPS.  It  tends  to  control  the  weight  assigned  to 
TOPS  relative  to  climatology  and  data. 

The  last  user-specified  parameter  used  in  OTIS  is 
the  noise-to-signal  ratio  for  the  subgrid  scale  error  Af 
defined  by 


L  ike  A,  B ,  and  C,  this  quantity  must  be  chosen  to 
reflect  the  statistics  of  ocean  thermal  variability.  For 
example.  A'  must  be  chosen  larger  in  dynamically 
active  regions,  where  subgrid-scale  eddies  produce 
large-amplitude  but  unresolved  perturbations  in  the 
thermal  field,  than  in  more  quiescent  regions  where 
the  eddy  field  is  weak.  Of  course,  A"  is  a  function 
of  the  grid  mesh  length,  becoming  smaller  as  the  grid 
becomes  finer  and  the  unresolved  portion  of  the 
thermal  field  decreases. 

Equation  (2-12)  is  used  along  with  the  RMS 
difference  between  bathy-observed  temperatures  and 
climatology  of  to  obtain  of  and  of.  Neglecting  tBATHY, 
which  is  small  compared  to  of  for  the  grids  currently 
utilized  by  OTIS,  (of)2  becomes 

(of)2  =  (of)2  +  (of)2  (2-13) 


or  just  the  sum  of  the  mean-square  resolvable  signal 
about  climatology  plus  the  mean-square  subgrid-scale 
noise.  From  (2-12)  and  (2-13), 


where  (ofv)2  is  the  mean-square  departure  of  the 
bathy  observations  from  climatology  and  Aff  is  the 
aveiage  value  of  Af,  both  calculated  over  the  entire 
analysis  domain.  The  data  used  in  this  calculation  con¬ 
sist  of  a  running  60-day  window  of  bathy  observations  # 

reported  to  FNOC  in  real  time,  amounting  to  about 
9000  reports. 

Using  the  assumption  of  equation  (2-16),  equations 
(2-14)  and  (2-15)  become 


(of)2  = 


K)2 

1  +  A". 


(2-17) 


(O2  =  TTr(°"‘):-  (2'18)  • 

With  Af  specified  by  the  user,  and  ohav  and  Aff 
calculated  internally  by  OTIS,  equations  (2-17)  and 
(2-18)  provide  of  and  of  for  use  in  equations  (2-4), 

(2-5),  (2-6),  (2-9),  and  (2-12).  Note  that  o*. ,  and  thus  f 

of  and  of,  are  determined  separately  for  each 
standard  level  of  the  OTIS  vertical  output  grid  (see 
Section  2.2). 

Note  also  from  equation  (2-17)  that  of  is 
independent  of  horizontal  location.  Thus,  implicit  in 
equation  (2-16)  is  the  assumption  that  variations  of  of  0 

with  location  are  produced  entirely  by  variations  in 

of- 

Equations  (2-2)  through  (2-12)  and  (2-17)  through 
(2-18)  constitute  a  closed  set,  which  can  be  solved  for 
the  weights  aki  and  ftk .  Once  these  weights  are 
obtained,  the  analyzed  temperature  is  calculated  from  # 

equation  (2-1).  Also,  from  the  basic  theory  of  OI  (see 
for  example  White  and  Bernstein,  1979),  the  RMS  error 
of  the  resolvable  thermal  field  of  can  then  be 
calculated  from 


(of)-  = 


K 


i  +  a: 


(2-14) 


(of)2  =  (of)2(l  -Pk-  I  **,%)• 


(2-19) 


(of)2  =  (2-15) 

For  practical  considerations  (e.g.,  bathy  data  cover¬ 
age),  it  is  not  possible  to  calculate  (of)2  directly  from 
bathy  data  at  every  observational  location.  Therefore, 
(of)2  is  estimated  from 


Recognizing  that  OTIS  initializes  the  TOPS  forecast, 

which  will  provide  Tk  for  the  following  day’s 

analysis,  the  value  of  (of)2  obtained  from  (2-19)  • 

becomes  (a k)jMT/AL  of  equation  (2-11)  (i.e.,  the  error 

in  the  TOPS  initial  conditions)  for  the  following  day’s 

analysis. 

Finally,  the  expected  departure  of  the  thermal  field 
from  an  independent  bathy  observation  made  precisely 
at  the  gridpoint  at  analysis  time  of  is  given  by  # 


„  ,  1  +  Af  .  , 

(o  )-  = - -  (o  ) 

'  I  +  A?  av 


(2-16)  of  =  [(of)2  +  (of)2  +  £2a4r„y]05  •  (2-20) 
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(2-23) 


The  three  terms  inside  the  brackets  represent  the 
contributions  of  analysis  error,  subgrid-scale  error, 
and  instrumental  error.  For  a  large  set  of  such 
comparisons,  approximately  68%  of  the  independent 
bathy  observations  would  differ  from  the  analyzed 
temperature  Tk  by  less  than  oek ,  and  approximately 
96%  would  differ  by  less  than  2 oek .  In  such 
dynamically  active  regions  as  the  Gulf  Stream  or  the 
Kuroshio,  ok  and,  hence,  ak  will  be  large  compared 
to  values  in  more  quiescent  regions,  such  as  the  eastern 
Pacific. 

2.2  Mixed-layer  Analysis 

Upon  conclusion  of  the  SST  analysis  at  gridpoint  k, 
OTIS  performs  an  OI  analysis  for  mixed-layer 
depth  (MLD).  In  analogy  to  equation  (2-1),  the 
expression  for  the  analyzed  MLD  K  « 

K  =  K  +  I  aki(h°  -  /*,f)  +  /w  -  V).  (2-21> 

I  =  1 

where  hk  is  the  climatological  MLD  at  the  gridpoint, 
hk  the  predicted  MLD  at  the  gridpoint  from  the 
24-hour  TOPS  forecast  from  the  previous  day’s  analysis, 
h"  the  observed  MLD  at  location  /,  h]  the 
climatological  MLD  at  location  i,  aki  the  weight 
applied  to  the  ith  observed  climatological  MLD 
anomaly  assimilated  at  the  gridpoint,  fik  the  weight 
applied  to  the  TOPS-predicted  MLD  anomaly  at  the 
gridpoint,  and  N  the  number  of  observations 
assimilated  at  the  gridpoint. 

The  climatological  MLD  hck  is  determined  from 
examination  of  the  climatological  temperature  profile, 
defined  at  gridpoint  k  on  the  vertical  grid  of  Table  2-1. 
Beginning  at  the  surface  and  working  downward, 
vertical  temperature  gradients  are  examined  to  find  the 
shallowest  pair  of  gridpoints  between  which  the  vertical 
gradient  is  less  than  -0.05°C  m  '.  The  layer  depth  hk 
is  then  taken  simply  as  the  depth  of  the  shallower  grid- 
point  of  this  pair,  and  is  thus  always  one  of  the  discrete 
levels  listed  in  Table  2-1.  Exactly  the  same  process  is 
carried  out  on  the  TOPS  and  bathy  profiles  to  define 
hk  and  h”.  The  TOPS  profile  is  also  defined  on  the 
upper  15  levels  of  Table  2-1,  but  the  depths  defining 
the  bathy  profiles  are  arbitrary;  thus,  h "  generally  falls 
somewhere  in  between  the  fixed  levels  of  Table  2-1. 

The  weights  akj  and  in  equation  (2-21)  are 
obtained  by  solving  the  N  +  1  equations  of  (2-2),  as 
before.  In  solving  (2-2),  the  parameters  r}lJ,  A",  and 
Xk  are  assumed  equal  to  the  corresponding  values  for 
temperature  analysis  at  the  100-m  depth.  This 
assumption  implies 

(2-22) 


(°k  )\ll  I)  ~  (K)(o'k  )|00m 
and 

(Omio  =  <*><o£>  .non,.  (2-24) 

where  the  subscript  “MLD”  indicates  values  are  for 
MLD  variations,  and  A-  is  a  constant. 

Although  hk  and  hk  can  only  take  on  the  discrete 
values  of  Table  2-1,  hk  resulting  from  the  solution 
of  (2-21)  and  (2-2)  can  take  on  any  value.  Thus,  it  can 
be  thought  of  as  a  “floating  level”  in  the  analysis 
system. 

With  the  analysis  for  hk  co\,,,.v.ie,  OTIS  sets  the 
temperatures  at  the  gridpomts  of  Table  2-1  that  are 
shallower  than  hk  according  to 

ttl 

(7?),,,  =  rn0  + 1  (Ar;)„ ,  (2-25) 

n  =  I 


Table  2-1.  Vertical  grid  for  OTIS  output. 


LEVEL 

DEPTH  IN  METERS 

1 

2.5 

2 

7.5 

3 

12.5 

4 

17.5 

5 

25.0 

6 

32.5 

7 

40.0 

8 

50.0 

9 

62.5 

10 

75.0 

11 

100.0 

12 

125.0 

13 

150.0 

14 

200.0 

15 

300.0 

16 

400.0 

17 

500.0 

18 

600.0 

19 

700.0 

20 

800.0 

21 

900.0 

22 

1000.0 

23 

1100.0 

24 

1200.0 

25 

1300.0 

26 

1400.0 

27 

1500.0 

28 

1750.0 

29 

2000.0 

30 

2500.0 

31 

3000.0 

32 

4000.0 

33 

5000.0 

(°>  Kit  d  -  (70(0,  )  ]00  m  , 


(2-28) 


where  ( Tk)m  is  the  analyzed  temperature  at  level  m, 
(r")0  is  the  analyzed  SST,  and  (A Tk)n  is  the  change 
in  TOPS-predicted  temperature  from  level  n  -  1  to  level 
n  (level  0  implies  the  surface).  Thus,  the  shape  of 
the  temperature  profile  from  the  surface  to  the  base 
of  the  mixed  layer  (e.g.,  isothermal,  weakly  stratified, 
or  multiple  thermocline)  is  controlled  exactly  by  TOPS. 

2.3  Sub-mixed-layer  Analysis 

After  the  mixed-layer  analysis  at  gridpoint  k  is 
complete,  OTIS  performs  an  analysis  for  the 
temperature  profile  below  the  analyzed  MLD.  Below 
this  level,  OTIS  utilizes  a  completely  variable  or 
“floating”  grid  system  designed  to  concentrate 
resolution  in  regions  of  high  vertical  gradient. 
Beginning  at  the  MLD  and  working  downward,  the 
mesh  spacing  of  this  floating  grid  A z  is  defined  as 

dT1'  d2Tp 

Az  =  50  m  +  (400m:C  ')— ^  -(lOm’C  2) - - 

oz  dz‘ 

(2-26) 

for  5  m  <  Az  <  80  m  , 

with  z  taken  positive  downward  from  the  sea  surface. 

With  the  vertical  grid  extending  from  the  base  of 
the  mixed  layer  downward  to  400  m,  defined  by  (2-26), 
OTIS  performs  an  analysis  for  temperature  at  each  of 
these  levels.  This  analysis  is  carried  out  by  solving  (2-1) 
and  (2-2)  at  each  of  the  floating  levels.  If  the  deep 
analysis  mode  is  selected,  OTIS  also  carries  out  the 
same  analysis  at  the  fixed  levels  17-33  of  Table  2-1. 

In  general,  the  bathy  observations  which  contribute 
to  the  sub-mixed-layer  analysis  extend  to  varying 
depths.  Bathy  observations  that  reach  to  less  than 
250  m,  or  which  have  ~  \(T"  -  T‘  ) |>  0.01  °C  m  1  near 
their  deepest  reported  depth  zH,  are  used  in  the 
analysis  only  for  levels  above  zB.  Observations 
extending  to  greater  than  250  m  depth  for  which 
ji_  7~‘‘ ) j  <  0.01  °Cm  1  near  zH  are  used  in  the 

analysis  at  all  levels  down  to  400  m,  however,  by  extrap¬ 
olating  (7"-  7~/  )  downward.  The  mean-square  error 
added  to  the  “observation”  below  zH  by  this  extrap¬ 
olation  is  estimated  from 

(o," ),;  v  =  (0.0225  +  0.02  exp  {  - 0.05  (zH  -250)}) 

(z  -  zB) , 

(2-27) 

where  z  is  the  depth  of  the  extrapolation  point  for  the 
/th  bathy  in  meters.  Equation  (2-7)  for  extrapolated 
bathy  observations  below  zH  then  becomes 


(°f  }  HA  IVY  ~  HATVY  +  (°/  )  +  (°,  )/:.V  • 

The  quantities  rj(y,  A",  and  A '’required  in  (2-2)  are 
derived  from  the  user-specified  parameters  Ak ,  Bk , 
Ck,  A  "and  fH1//n  as  before.  In  calculating  A  "and 
A£  via  (2-4),  (2-5),  and  (2-7),  the  quantities  o‘,  o'1 
and  of  are  vertically  interpolated  from  the  fixed- 
levels  of  Table  2-1  to  the  floating  levels.  Equation 
(2-28)  is  used  in  place  of  (2-7)  for  extrapolated  bathy 
observations  below  zH  as  discussed  above.  Note  that 
the  parameters  Ak ,  Bk,  Ck,  and  A"  for  sub-mixed- 
layer  analysis  can  be  assigned  values  by  the  user 
distinctly  different  from  those  at  the  surface.  This  is 
particularly  important  for  Ck ,  which  generally  has  a 
much  smaller  value  below  the  MLD  reflecting  the  much 
longer  time-scale  variability  there. 

Upon  conclusion  of  the  sub-mixed-layer  temperature 
analysis,  the  analyzed  temperatures  at  the  floating 
levels  defined  by  (2-26)  are  vertically  interpolated  to 
the  fixed  levels  of  Table  2-1  between  the  MLD  and  the 
400-m  depth  for  output.  Thus,  the  final  OTIS  surface- 
to-400  m  (shallow  mode)  or  surface-to-5000  m  (deep 
mode)  temperature  profile  is  defined  at  the  levels  of 
Table  2-1.  The  upper  15  levels  of  this  grid  are  shared 
with  TOPS,  while  the  low'er  18  are  shared  with  the 
GDEM  ocean  thermal  climatology  (see  Section  8). 

As  a  final  output  step  for  the  deep  analysis  mode, 
the  analyzed  vertical  temperature  profile  below  400  m 
is  altered  where  necessary  to  blend  with  the  vertical 
temperature  profile  analyzed  above  using  the  Navy- 
standard  merge  algorithm.  This  step  is  necessary 
because  400  m  represents  a  discontinuity  in  data 
density,  corresponding  to  the  nominal  depth  of  most 
bathy  observations. 


3.0  User-specified  Parameters 

The  nine  user-specifed  parameters  used  in  OTIS  are 
listed  in  Table  3-1 . 


Table  3-1.  User-specified  parameters. 


PARAMETER 

VALUE 

£SHIP 

2.29°C 

£  BATHY 

0.2CTC 

£mcsst 

0.71  "C 

eMCSST 

0.2  2”C 

C 

1/(20  days)  above  the  MLD.  1/(60  days)  below 
the  MLD. 

A 

Variable,  defined  in  latitude/longitude 
bands,  separate  values  for  above  and 
below  the  MLD 

8 

G 

A" 

The  values  of  £s///p  and  £BATHY  are  obtained  from 
Earle  (1985)  and  from  White  and  Bernstein  (1979), 
respectively.  The  values  of  tMCSST  and  eKfCSST  are 
based  on  Hawkins  et  al.  (1986)  and  Strong  and 
McClain  (1985).  The  values  of  A  and  B  range  from 
1/(2500  km)  in  low-latitude  open-ocean  regions  to 
1/(100  km)  in  coastal  regions.  The  values  of  G  and  A" 
vary  widely.  In  general,  G  is  smaller  in  midlatitude 
open-ocean  regions  than  in  tropical  or  western¬ 
boundary  current  regions,  where  the  atmospheric 
forcing  and/or  physics  of  TOPS  is  not  fully  adequate. 
In  general,  A”  is  smaller  in  quiescent  regions  (e.g., 
eastern  Pacific)  than  in  dynamically  active  regions 
characterized  by  strong  mesoscale  eddies  which  are 
unresolved  by  the  horizontal  grid  (e.g.,  Kuroshio  or 
Gulf  Stream  regions). 

4.0  Data  Screening 

OTIS  utilizes  only  data  that  fall  within  certain 
temporal  and  spatial  windows  (see  Phoebus,  1988).  For 
the  SST  analysis,  the  temporal  window  includes  all  data 
0  to  60  hours  prior  to  the  OTIS  run  time.  For  the  MLD 
and  sub-mixed-layer  analyses,  the  window  includes  all 
bathys  from  0  to  60  days  prior  to  the  OTIS  run  time. 

The  spatial  window  is  defined  as  ±  1  correlation  scale 
in  the  north-south  and  east-west  directions  about  the 
gridpoint.  Before  data  are  selected  from  this  window, 
however,  the  bathy  data  are  screened  by  applying  a 
sorting  process  that  limits  the  bathy  data  density  in 
2.5°  x  2.5°  latitude/longitude  squares  to  the  60  most 
recent  reports,  with  no  more  than  30  taken  from  the 
same  platform.  This  thinning  process  is  a  practical 
necessity  for  handling  ocean  weather  stations,  which 
typically  report  two  bathy  observations  per  day. 

In  addition,  OTIS  applies  a  ship-tracking  algorithm, 
patterned  after  the  approach  of  Mendenhall  and 
Cuming  (1980),  to  screen  out  bathy  reports  with 
obvious  position  errors.  These  types  of  errors  are  fairly 
common  in  the  real-time  bathy  data  set  and  result 
primarily  from  coding  and  radio  transmission  errors 
(e.g.,  latitude  and  longitude  transposed,  wrong 
hemisphere,  etc.).  Thus,  OTIS  monitors  the  successive 
positions  of  each  ship  reporting  bathy  observations  at 
any  time  during  the  previous  60-day  period  and  flags 
movements  that  reflect  impossibly  large  ship  speeds. 
Bathy  observations  corresponding  to  the  erroneous 
position  reports  are  excluded  from  further  processing 
by  OTIS. 

For  the  SST  analysis,  up  to  700  reports  within  the 
spatial  window  are  collected  without  regard  to  relative 
correlation  or  data  type  (i.e.,  ship,  bathy,  MCSST). 
Data  in  excess  of  the  700-report  limit  are  ignored,  but 
this  rarely  happens  in  practice.  These  data  are  then 
sorted  to  select  the  15  reports  most  highly  correlated 
with  the  gridpoint,  again  without  regard  to  data  type. 


The  15  most  highly  correlated  reports  are  then 
screened  for  spurious  data  using  the  “buddy  check” 
approach  of  Phoebus  (1989),  which  is  patterned  after 
that  of  Dimego  et  al.  (1985).  Hence,  the  inequality 

\iK-  t; )  -  (T“  -  r/)|  <  <3  - 1.5  %)o'  (4-i) 

is  examined.  If  the  inequality  is  satisfied,  then 
observations  /  and  j  corroborate  each  other  and  are 
thus  both  given  “keep  flags”  equal  to  the  correlation 
between  them,  r)u;  if  the  inequality  is  not  satisfied, 
then  the  observations  contradict  each  other  and  are  thus 
both  given  “toss  flags”  equal  to  r\j  .  Following 
examination  of  all  remaining  pairs  of  reports,  the  single 
observation  with  the  highest  summation  of  toss  flags 
is  removed  if  the  summation  exceeds  2.  Observations 
whose  keep  flags  sum  to  more  than  2  are  retained 
regardless  of  toss  flag  settings.  This  process,  beginning 
with  the  examination  of  the  inequality  (4-1),  is  then 
repeated  until  no  observations  have  toss  flags  that  sum 
to  more  than  2.  The  remaining  observations  are  utilized 
in  the  SST  analysis  at  the  gridpoint. 

Screening  of  the  bathy  data  used  in  the  MLD  and 
sub-mixed-layer  analyses  proceeds  in  similar  fashion. 
Up  to  700  bathy  reports  are  selected  in  the  spatial 
window  without  regard  to  their  relative  correlation  with 
the  gridpoint.  The  data  are  then  sorted  to  retain  only 
the  15  observations  most  highly  correlated  with  the 
gridpoint.  Finally,  the  buddy  check  of  equation  (4-1) 
is  applied  to  the  15  most  highly  correlated  observations 
to  eliminate  spurious  reports. 

OTIS  also  compiles  a  list  of  bathy  reports  rejected 
via  the  buddy  check  described  above.  This  list  can  be 
sorted  in  various  ways  and  output  daily  or  weekly  as 
an  aid  to  analysts  performing  subjective  bathy  error 
correction  via  the  Quality  Improvement  Profile  System 
(QUIPS)  at  FNOC  (see  McLain  et  al.,  1986). 

OTIS  also  includes  the  capability  to  exclude 
individual  bathy  reports  identified  by  quality  control 
personnel  as  being  suspect,  even  if  they  pass  the  internal 
OTIS  data  screening  checks.  Information  identifying 
these  reports  is  simply  typed  into  a  disk  file,  which 
OTIS  automatically  checks  each  run.  This  process 
allows  convenient  application  of  an  additional  level  of 
data  quality  control. 

Finally,  an  additional  data  screening  capability  can 
be  used  in  certain  regional  applications  of  OTIS.  An 
option  exists  to  decorrelate  data  from  the  analysis 
across  fronts  and  eddy  boundaries,  which  are 
determined  from  an  analyst’s  interpretation  of  satellite 
imagery  and  represented  on  the  FNOC  data  base  in 
message  format.  If  this  option  is  selected,  then  OTIS 
will  read  the  front  and  eddy  messages  from  the  FNOC 
data  base,  determine  the  location  of  fronts  and  eddy 
boundaries,  and  truncate  the  spatial  windows  used  for 


data  selection  around  each  gridpoint  to  ensure  that  data 
on  one  side  of  a  front  or  eddy  boundary  does  not 
contribute  to  the  analysis  at  gridpoints  on  the  other 
side  of  the  front  or  eddy  boundary.  Thus,  in  effect, 
a  water-mass-based  analysis  is  achieved  with 
observations  made  in  a  particular  water  mass  (e.g., 
shelf  water,  slope  water,  Gulf  Stream  water.  Sargasso 
Sea  water,  warm-core  ring,  cold-core  ring,  etc.) 
affecting  only  gridpoints  contained  within  that  water 
mass. 


5.0  Prediction  Models 

The  predicted  temperature  T%  required  by  OTIS  is 
normally  taken  from  a  24-hour  forecast  produced  by 
the  TOPS  mixed-layer  model  (see  Clancy  and  Martin, 
1981;  Clancy  and  Poliak,  1983;  Martin  et  al.,  1985) 
initialized  from  the  previous  day’s  OTIS  analy¬ 
sis.  Thus,  OTIS  and  TOPS  are  coupled  in  an 
analysis-forecast-analysis  data  assimilation  cycle.  OTIS 
also,  however,  includes  the  capability  to  use  a 
“climatological  trend”  model  in  place  of  TOPS  as 
described  below. 

5.1  TOPS  Model 

The  aim  of  TOPS  is  to  forecast  the  response  of  the 
mixed  layer  and  mid-ocean  fronts  to  atmospheric 
forcing.  The  prognostic  equations  for  TOPS  express 
conservation  of  temperature,  salinity,  and  momentum 
in  the  upper  ocean  as  follows: 


dT 

dt 


.|(,n.|KntA(0  +  0) 


1  dF  d 


(5-1) 


where  7" is  the  temperature,  S  the  salinity,  u  and  v  the 
x-  and  y-components  of  the  current  velocity  (\  and  y 
horizontal  coordinates  defined  relative  to  the  TOPS 
grid),  w  the  z-component  of  the  current  velocity,  /-'the 
downward  flux  of  solar  radiation,  p  a  reference  density, 
c  the  specific  heat  of  seawater,  d  a  damping  coefficient, 
x  a  diffusion  coefficient, /the  Coriolis  parameter,  A 
the  horizontal  eddy  diffusion  coefficient,  /  the  time, 
z  the  vertical  coordinate  (positive  upward  from  the  sea 
surface),  and  ua,  vu,  wa  the  x-,  y-,  and  z-components 
of  an  advection  current  defined  below. 

The  terms  involving  d  in  (5-3)  and  (5-4)  represent 
the  drag  force  caused  by  the  radiation  stress  at  the  base 
of  the  mixed  layer  associated  with  the  propagation  of 
internal  wave  energy  downward  and  away  from  the 
wind-forced  region  (see  Pollard  and  Millard,  1970). 
The  terms  involving  x  account  for  very  weak 
“background”  eddy  diffusion  (due  to  intermittent 
breaking  of  internal  waves,  for  example)  that  exists 
even  below  the  mixed  layer.  The  terms  involving  A 
represent  horizontal  eddy  diffusion  patterned  after  that 
of  Haney  (1974). _ _ 

The  terms  w'7','  w'S,'  w'u,'  and  w’ v'  represent 
the  vertical  eddy  (i.e.,  turbulent)  fluxes  of  temperature, 
salinity,  and  momentum.  These  terms  account  for  the 
major  portion  of  the  mixed-layer’s  response  to 
atmospheric  forcing  and  are  parameterized  with  the 
Level-2  turbulence  closure  theory  of  Mellor  and 
Yamada  (1974)  and  Mellor  and  Durbin  (1975).  The 
upper  boundary  conditions  for  these  terms  constitute 
the  surface  fluxes  of  heat,  moisture,  and  momentum, 
which  drive  the  TOPS  forecasts.  Thus, 


dT , 


[-w'T  +  x  -^-1 

az  z  =  o 
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as. 


■w'S'  +  x-^J 

OZ  :  =  0 


du. 


QoSo 
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l-w'u'  +  X 

dz  z  =  o 
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and 
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where  Ba  is  the  surface  infrared  radiative  heat  flux, 
Ho  is  the  surface  sensible  heat  flux,  LQa  is  the  surface 
latent  heat  flux,  Fa  is  the  surface  solar  heat  flux,  SQ 
is  the  surface  salinity,  pa  is  the  density  of  the  air  near 
the  surface,  CD  is  a  drag  coefficient,  and  U  and  V  are 
the  components  of  the  wind  velocity  vector  at  a 
reference  height  of  19.5  m  above  the  sea  surface.  The 
quantities  Bu,  Ha ,  LQ0,  Fa,  U,  and  V  are  provided 
by  either  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS;  see  Rosmond,  1981)  or 
the  Navy  Operational  Regional  Atmospheric  Prediction 
System  (NORAPS;  see  Hodur,  1982)  models, 
depending  on  the  scale  and  area  of  application. 

Because  there  are  no  horizontal  pressure  gradient 
terms  in  (5-3)  and  (5-4),  u  and  v  represent  the  wind- 
drift  (i.e. ,  Ekman  plus  inertial)  component  of  the 
current.  Neglect  of  horizontal  pressure  gradients  here 
is  motivated  by  the  fact  that  geostrophic  currents 
generally  do  not  play  an  important  role  in  mixed-layer 
dynamics.  Still,  TOPS  allows  geostrophic  currents 
supplied  from  an  external  source  to  contribute  to  the 
advection.  Thus, 


n'  =  (T“k),smAl  +  A  Tcim  (5-13) 

where  (Tk),smAL  is  the  temperature  at  the  appropri¬ 
ate  level  resulting  from  the  previous  OTIS  analysis 
(typically  24  hours  earlier)  and  A Tcint  is  the 
climatological  change  in  temperature  at  this  level 
between  the  time  of  the  previous  analysis  and  the 
current  analysis  time.  As  is  the  case  for  TOPS,  the  error 
growth  rate  Gk  is  specified  in  latitude-longitude 
rectangles  by  the  user. 

This  option  provides  the  flexibility  to  run  OTIS  in 
regions  for  which  there  is  no  TOPS  model  available. 

6.0  Idealized  Cases 

In  this  section  we  present  idealized  cases  that 
illustrate  certain  aspects  of  the  OI  analysis  technique 
as  implemented  in  OTIS. 

6.1  Data-void  Gridpoint 

In  the  case  of  a  data-void  gridpoint,  equation  (1) 
reduces  to 


(5-10)  T*  =  Tc  +  pk(T[-  Tk),  (6-1) 

(5-H)  which  can  be  written 


and 


Tak  =  7/  -  vk(Tkp  -  7p  , 


+ 


a 

dyVa 


)dz , 


(5-12) 


where 

vk=  1  ~Pk- 


(6-2) 


(6-3) 


where  ug  and  vg  are  the  x-  and  y-components  of  the 
geostrophic  current.  The  source  of  these  geostrophic 
currents  will  ultimately  become  eddy-resolving 
dynamical  ocean  circulation  models  (collectively 
identified  by  the  acronym  OCEANS)  that  will  run 
independently  of  TOPS  (see  Hurlburt,  1984).  This 
input  will  allow  a  much  improved  representation  of 
the  thermal  field  associated  with  fronts  and  eddies,  even 
in  the  absence  of  in  situ  data.  Even  without  geostrophic 
advection  currents,  however,  TOPS  contains  the 
essential  physics  necessary  to  predict  much  of  the 
evolution  of  the  subarctic  and  subtropical  fronts 
(Roden,  1977;  Roden  and  Paskausky,  1978). 

See  Clancy  and  Poliak  (1983)  and  FNOC  (1986)  for 
a  detailed  and  complete  discussion  of  TOPS  and  its 
implementation  at  FNOC. 


The  second  term  on  the  right  side  of  (6-2)  is  an 
“anomaly  decay  term,”  which  attempts  to  return  the 
analyzed  temperature  to  climatology.  If  the  TOPS 
predicted  temperature  Tk  does  not  change  from  day 
to  day,  then  the  analyzed  temperature  Tk  will  decay 
toward  climatology  Tk,  with  vk  representing  the 
percent  anomaly  reduction  per  analysis  cycle  (i.e.,  per 
day).  The  number  of  days,  D,  required  to  reduce  the 
anomaly  to  10%  of  the  starting  value  is  thus 


D  = 


In  (0.1)  =  In  (0.1) 

ln(\-vky  lln((ik ) 


(6-4) 


Equation  (2-2)  reduces  to 

1=(1+  Xp)ftk  .  (6-5) 


5.2.  Climatological  Trend  Model 

OTIS  has  the  capability  to  utilize  a  “climatological 
trend”  model  in  lieu  of  TOPS  to  provide  Tp  at  any 
or  all  of  the  OTIS  analysis  levels  listed  in  Table  2-1. 
In  this  case, 


From  (6-5)  and  (2-5),  ftk  can  be  expressed  in  terms  of 
ock  and  opk.  Thus, 


Pk  = 


(ock)2 


«)2  +  (o/V 


(6-6) 
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From  (6-6)  and  (2-19),  the  mean-square  error  of  the 
resulting  analysis  can  be  written  as 

,  (0*0/') 2 

(op-  =  -2  •  (6-7) 

(0*)'  +  K> 


no  restoration  toward  climatology.  If  Gk  equals  olk, 
then  v*  equals  approximately  0.62  and  (6-2) 
reduces  to 


Tak  =  Tpk  -  0.62  (7?  -  Tf)  . 


(6-12) 


Now  ,  with  no  data  inputs,  the  mean-square  analysis 
error  will  reach  a  steady  state  (op2  at  gridpoint  k. 
But,  from  (2-11),  recognizing  that  (ok)jMTJAL  will  be 
given  by  (op\  we  have 


In  this  case,  a  rather  strong  restoring  force  toward 
climatology  is  imposed  on  each  analysis  cycle.  As  Gk 
approaches  infinity,  vk  approaches  1  and  (6-2) 
reduces  to 


(0p‘  +  g:  . 


'J'Q  _  ’T'C 

lk  ~  1 k  ‘ 


(6-13) 


Substitution  of  (6-8)  into  (6-7)  and  replacement  of 
(ok ) *  in  (6-7)  with  ( ok )2,  yields  a  quadratic  equation 
for  (op2,  whose  solution  is 


(0*')2  =  Gp^G2  +  K)2]05  -  \  G2k 


Certain  limiting  values  for  the  TOPS  error  growth 
parameter  are  of  interest  in  (6-9).  If  Gk  equals  0,  then 
(op2  also  equals  zero.  Thus,  if  the  prediction  model 
is  perfect  (i.e.,  has  zero  error  growth  rate),  then  the 
analysis  error  tends  toward  zero  in  the  asymptotic 
steady  state.  If  Gk  equals  ak,  then  (op2  equals 
approximately  0.62 (op2.  Thus,  even  with  this  large 
error  growth  rate,  the  prediction  model  is  able  to 
contribute  some  useful  information  to  the  analysis; 
w  ithout  its  input,  the  analysis  error  would  tend  toward 
(op2  in  the  steady  state.  If  Gk  equals  3oP  then  (op2 
equals  approximately  0.91  (op2,  and  the  TOPS 
predictions  contribute  little  to  reducing  the  analysis 
error.  Finally,  as  Gk  approaches  infinity,  (op2 
approaches  (op2. 

From  (6-3),  (6-6),  (6-8),  and  (6-9),  the  anomaly 
reduction  factor  vk  can  be  expressed  as 

ic2  +  GjiG2  +  (op2]0-5 

(Op2  +  \  Gl  +  epic2  +  (op2]0-5 

(6-10) 


Limiting  values  for  Gk  are  again  of  interest.  If  Gk 
equals  0,  then  vk  also  equals  0  and  (6-2)  reduces  to 


7-  W  __ 

1  k  ~  1  k 


(6-11) 


Thus,  if  the  prediction  model  has  no  skill  at  all 
(Gk  —  °°),  the  analysis  simply  reproduces  climatology 
in  the  long-term  absence  of  data,  regardless  of  what 
changes  the  prediction  model  forecasts.  Table  6-1 
shows  vk  and  D  as  a  function  of  Gk  for  ak  equal  to 
1.0°C. 

From  (2-19)  and  (6-3),  the  steady-state  mean-square 
error  of  the  analysis  can  be  expressed  as 


(op2  =  (op2(l  -f}k)  =  ( op2vA.  . 


(6-14) 


Table  6-1.  Anomaly  reduction  factor  vk  and  number  of  days  D  for 
the  analyzed  thermal  anomaly  to  be  reduced  to  10%  of  initial  value 
at  data-void  gridpoint  as  a  function  of  TOPS  error  growth  rate  Gk . 
Assumes  one  analysis  cycle  per  day  and  ok  equal  to  1.0°C. 


Thus,  in  the  case  of  a  perfect  model,  the  analysis  will 
simply  carry  the  model  prediction  forward  in  time  with 


Gk(°  C) 

vk 

D  (days) 

0.00 

0.00 

00 

0.03 

0.03 

76.7 

0.04 

0.04 

57.6 

0.05 

0.05 

46.0 

0.06 

0.06 

38.4 

0.07 

0.07 

33.0 

0.08 

0.08 

28.8 

0.09 

0.09 

25.6 

0.10 

0.09 

23.0 

0.15 

0.14 

15.4 

0,20 

0.18 

11.5 

0.25 

0.22 

9.2 

0.30 

0.26 

7.7 

0.35 

0.29 

6.6 

0.40 

0.33 

5.8 

0.60 

0.45 

3.9 

1.00 

0.62 

2.4 

2.00 

0.83 

1.3 

3.00 

0.91 

1.0 

00 

1.00 

0.0 

and 


I 

Thus,  the  values  of  vk  in  Table  6-1  can  also  be 
interpreted  as  the  ratio  of  steady-state  mean-square 
analysis  error  to  the  mean-square  error  of  climatology 
for  various  TOPS  error  growth  rates  Gk . 

*  6.2  Isolated  Observation  at  a  Gridpoint 

Here  we  investigate  the  response  of  the  analysis 
system  to  an  isolated  SST  observation  following  an 
extended  period  for  which  no  data  were  available.  In 
this  situation,  equation  (2-1)  becomes 

1 

n  =  rk  +  o,(r;  -  T‘  )  +  (lk  (Tk  -  T‘  ),  (6-15) 


ftk  = 


>1 " 
aa- 


K 


+  xpx‘‘  +  a; 


(6-22) 


Four  cases  are  of  interest. 

Case  1:  Perfect  prediction,  imperfect  observation. 
As  demonstrated  in  the  previous  section,  a  perfect 
prediction  model  (i.e. ,  one  for  which  Gk  =  0)  yields 
oak  =  osk  =  0,  thus  implying  Xpk  =  0.  From  (6-15), 
(6-21),  and  (6-22)  with  Xpk  =  0  and  X°k  #  0,  we  find 


which  can  be  written 


'pa 

1  k 


' pP 

1  k  * 


(6-23) 


Tk  =  <*JZ  +  ft  Jr  +  yJr  (6-16) 

where  yk  is  the  weight  applied  to  climatology, 
given  by 

Yk=  1  -  -  ftk  •  (6-17) 

Furthermore,  equation  (2-2)  reduces  to 


»J*,  =  0  +  A  °k)ak  +  rikiftk,  (6- 18a) 

and 


1  -  H*,0*  +  (1  +  * k)ftk  > 


(6- 18b) 


where  we  have  made  use  of  the  fact  that  rjik  =  r]kr 
Simultaneous  solution  of  (6- 18a)  and  (6- 18b)  yields 


»?*,■** 


(i  -  ik,)  +  K  +  KK  +  K 


(6-19) 


and 


ftk 


i 


_ (i  -  vjkft  +  K _ 

(1  -  r)2ki)  +  \p  +  \p\%  +  \° 


(6-20) 


6.2.1  Initial  Response  of  the  Analysis  to  an  Observation 
Made  Precisely  at  the  Gridpoint  at  Analysis  Time 

If  the  isolated  observation  is  made  precisely  at  the 
gridpoint  at  analysis  time,  then  r\kj  =  1  and  (6-19) 
I  and  (6-20)  reduce  to 


= 


K 


\p  +  xpx°k  +  x°k 


(6-21) 


Thus,  the  data  are  completely  ignored  and  the  analysis 
reproduces  the  TOPS  prediction  exactly. 

Case  2:  Imperfect  prediction,  perfect  observation. 

In  this  case,  the  observation  is  assumed  perfect;  and 
thus,  A£  =  0.  Note  that  this  case  could  occur  only  if 
both  the  instrumental  and  the  subgrid-scale  errors  of 
the  observation  were  zero.  From  (6-15),  (6-21),  and 
(6-22)  with  A"  =  0  and  A£  =£  0,  we  find 

T°k  =  Tk  .  (6-24) 

Thus,  the  TOPS  prediction  is  ignored  completely,  and 
the  analysis  system  reproduces  the  observation  exactly. 

Case  3:  Perfect  prediction,  perfect  observation. 

In  this  case,  Xpk  =  X°k  =  0.  Also,  since  the 
prediction  and  the  observation  are  both  “perfect,”  they 
must  equal  each  other.  Thus,  Tk  =  Tk.  From  (6-15) 
with  Tk  =  Tk,  and  (6-18a)  or  (6-18b)  with 
A*  =  A£  =  0  and  =  1,  we  find 

Tak  =  Tk  =  Tp  .  (6-25) 

Thus,  the  analysis  exactly  reproduces  the  observation 
and  the  TOPS  prediction,  which  is  equal  to  the 
observation. 

Case  4:  Imperfect  prediction,  imperfect  observation. 

This  case  is,  of  course,  the  general  “real-world”  case 
for  which  Xk  and  A£  both  differ  from  zero.  Equations 
(2-4)  through  (2-9),  (2-12),  (6-8),  (6-9),  (6-21),  (6-22), 
and  (6-17),  along  with  the  instrumental  error 
parameters  of  Table  3-1,  yield  the  weight  applied  to 
the  isolated  observation  ak ,  the  weight  applied  to  the 
TOPS  prediction  pk,  and  the  weight  applied  to 
climatology  yk  as  a  function  of  Gk,  Xnk,  olk,  and  the 
observational  type  (an  additional  parameter,  M,  the 
number  of  MCSST  reports  averaged  to  form  a  super 
ob,  is  also  a  factor  in  determining  the  weight  for  the 
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MCSST  data).  Tables  6-2  and  6-3  show  ak . ,  ftk,  and 
yk  for  the  different  observational  types  as  a  function 
of  Gk  and  A"  for  ok  =  1.0  °C. 

As  illustrated  by  the  data  presented  in  Tables  6-2 
and  6-3,  the  weight  given  to  the  TOPS  predictions  ftk 
falls  off  rapidly  with  increasing  TOPS  error  growth 
rate  Gk.  The  analysis  effectively  “loses  confidence” 
in  the  TOPS  predictions  and  relies  increasingly  on  the 
data  and  climatology  (i.e.,  ak  and  yk  increase)  as  the 
TOPS  error  growth  rate  increases.  Increasing  ynk 
implies  less  accurate  data  as  a  result  of  higher  subgrid- 
scale  noise  contamination  in  the  observations.  Thus, 
for  a  fixed  error  growth  rate,  the  weight  given  to  TOPS 
and  the  weight  given  to  climatology  yk  increase  with 
increasing  noise-to-signal  ratio  A"  as  the  weights  given 
to  the  various  types  of  observations  ak  decrease.  For 
M  (the  number  of  MCSST  reports  averaged  to  form 
a  super  ob)  equal  to  1,  the  highest  weight  is  assigned 
to  bathy  reports,  the  next  highest  weight  to  MCSST 
reports,  and  the  least  weight  to  ship  reports,  reflecting 
the  relative  instrumental  accuracies  of  these  data.  As 
A"  increases  for  a  fixed  value  of  Gk,  the  weights 


assigned  to  the  observations  tend  to  become  more  even 
as  the  subgrid-scale  noise  obscures  the  differences  in 
instrumental  accuracy.  For  M  =  50  and  M  =  100,  the 
MCSST  super  obs  become  the  dominant  observation. 
Note  that  increasing  M  from  50  to  100  has  only  a  small 
impact  on  the  weight  assigned  to  the  MCSST  super  ob. 

6.2.2  Time-varying  Response  of  the  Analysis 
to  an  Observation  Made  Precisely  at  the 
Gridpoint  at  Analysis  Time 

In  this  section  we  investigate  the  time-varying 
response  of  the  analysis  system  to  an  isolated  bathy 
SST  observation  made  precisely  at  the  analysis  grid- 
point.  Using  equations  (2-3),  (2-4),  (2-5),  (2-7),  (2-11), 
(2-12),  (2-19),  (6-9),  (6-14),  (6-17),  (6-19),  and  (6-20), 
we  produce  time  series  of  ak,  ak,  (3k,  and  yk,  which 
illustrate  this  response. 

Figure  6- 1(a)  shows  time  series  of  the  weight  assigned 
to  the  observation  crk,  the  weight  assigned  to  TOPS 
1 3k ,  and  the  weight  assigned  to  climatology  yk  for 
A£  =  0.5,  ok  =  1°C,  and  Gk .  =  0.2°C.  On  Days  1 


Table  6-2  Weights  given  to  isolated  observation  at  a  gridpoint  ok,  TOPS  prediction  ftk,  and  climatology  yk  as  a  function  of  TOPS  error 
growth  rate  G„  and  observational  type,  assuming  kk  =  0.5  and  ok  =  1.0°C.  M  is  the  number  of  MCSST  reports  averaged  at  a  gridpoint  to 
form  a  super  ob 


G*(°C) 

a,  ft,  y  for  Ship  Ob 

a.  ft,  y  for  Bathy  Ob 

a.  p.  y  for  MCSST  Super  Ob 

0.10 

0.02  0.89  0.09 

0.1 6  0.76  0.08 

0.05 

0.86 

0.09 

(M 

=  1) 

0.52 

0.43 

0.05 

(M 

=  50) 

0.59 

0.37 

0.04 

(M 

=  100) 

0.20 

003  0.79  0.18 

0.25  0.61  0.14 

0.09 

0.75 

0.16 

(M 

=  1) 

0.67 

0.27 

0.06 

(M 

=  50) 

0.72 

0.23 

0.05 

(M 

=  100) 

1.00 

0.10  0.34  0.56 

0.53  0.18  0.29 

0.25 

0.29 

0.46 

(M 

=  1) 

0.87 

0.05 

0.08 

(M 

=  50) 

0.90 

0.04 

0.06 

(M 

=  100) 

3  00 

0  14  0.08  0  78 

0.63  0.03  0.34 

0.32 

0.06 

0.62 

(M 

=  1) 

091 

0.01 

0.08 

(M 

=  50) 

0.93 

0.01 

0.06 

(M 

=  100) 

Table  6-3.  Same  as  Table  6-2,  but  for  kk  =  1-5. 


G*(°C) 

a.  p.  y  for  Ship  Ob 

a.  P.  y  for  Bathy  Ob 

a,  p,  y  for  MCSST  Super  Ob 

0. 1C 

0.01  0.89  0.10 

0.06  0.85  0.09 

0.03 

0.87 

0.10 

(M 

=  1) 

0.47 

0.47 

0.06 

(M 

=  50) 

0.55 

0.40 

0.05 

(M 

=  100) 

0  20 

0.03  080  0.18 

0.11  0,73  0.16 

0.06 

0.77 

0.17 

(M 

=  1) 

0.62 

0.31 

0.07 

(M 

=  50) 

0.69 

0.25 

0.06 

(M 

=  100) 

1  00 

0.08  0.35  0.57 

0.29  0.27  0.44 

0.18 

0.31 

0.51 

(M 

=  1) 

0  85 

0.06 

0.09 

(M 

=  50) 

0.89 

0.04 

0.07 

(M 

=  100) 

3  00 

0.12  0.08  0.80 

0.37  0.06  0.57 

0.24 

0.07 

0.69 

(M 

=  1) 

0.89 

0.01 

0.10 

(M 

=  50) 

0.92 

0.01 

0.07 

(M 

=  100) 
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Figure  6 -1(a).  Time  series  of  the  weight  given  to  TOPS  ftk ,  the 
weight  given  to  the  observation  ok,  and  the  weight  given  to 
climatology  yk  for  a  single  bathy  SST  observation.  The 
observation  was  made  precisely  at  the  analysis  gridpoint  at 
Day  3  following  an  extended  period  for  which  no  data  were 
available.  The  noise-to-signal  ratio  Ak,  the  RMS  variation 
of  the  resolvable  thermal  anomalies  about  climatology  ok ,  and 
the  TOPS  error  growth  rate  Gk  are  taken  to  be  0.5,  I°C, 
and  0.2°C,  respectively. 


DAY 


Figure  6-Kb).  Same  as  Figure  6-l(a),  but  for  Gk  =  I.O°C. 


and  2,  ftk  and  yk  exhibit  their  data-void  steady-state 
values  associated  with  the  prescribed  value  of  Gk.  The 
bathy  observation  is  made  on  Day  3  and  the  weights 
immediately  adjust  in  response  to  it.  The  magnitude 
of  pk  drops  abruptly  in  response  to  the  observation, 


and  then  increases  asymptotically  back  to  its  data-void 
steady-state  value.  This  change  reflects  the  updating 
of  the  TOPS  prediction  with  the  information  supplied 
by  the  observation  at  Day  3,  and  the  carrying  of  this 
information  forward  in  lime  by  the  TOPS  model.  The 
magnitude  of  ak  decays  monotonically  to  zero  fium 
its  initial  non-zero  value  at  Day  3,  reflecting  the  loss 
of  information  as  the  observation  “ages.”  The  weight 
given  to  climatology  yk  takes  on  a  minimum  at  about 
Day  13  before  returning  asymptotically  to  its  data-void 
steady-state  value. 

Figure  6- 1(b)  shows  time  series  of  cr*.,  ftk,  and  yk  for 
the  same  scenario  as  Figure  6- 1(a),  except  with 
Gk  =  1.0°C.  In  this  case,  the  analysis  has  much  less 
“confidence”  in  the  TOPS  prediction  because  of  the 
larger  TOPS  error  growth  rate.  As  a  result,  compared 
to  the  previous  case,  more  weight  is  given  to  both  the 
observation  and  climatology,  and  less  weight  is  given 
to  TOPS.  Also,  the  weights  return  to  their  data-void 
steady-state  values  more  quickly  in  this  situation  than 
before,  as  TOPS  is  less  effective  in  carrying  infor¬ 
mation  forward  in  time. 

6.3  Two  Observations  Near  a  Gridpoint 

In  this  section  we  investigate  the  initial  response  of 
the  analysis  system  to  two  bathy  SST  observations 
made  precisely  at  analysis  time  following  an  extended 
period  for  which  no  data  were  available.  To  simplify 
the  geometry,  the  observations  are  distributed  in  a 
straight  line  that  passes  through  the  gridpoint  location. 
By  varying  the  location  of  the  observations  along  this 
line,  we  can  illustrate  the  dependence  of  the  assigned 
weights  and  the  resulting  analysis  error  on  data 
distribution. 

The  weights  given  to  the  two  observations,  a,  and 
a-,,  and  the  weight  given  to  TOPS  ft  are  determined 
from  (2-2),  (2-3),  (2-4),  (2-5),  (2-7),  (2-11),  (2-12),  and 
(6-9).  The  analysis  error  oak  is  determined  from  (2-19). 
In  the  cases  that  follow,  we  assume  \nk  =  0.5, 
Gk  =  0.2°C,  ok  =  1°C,  and  Ak  =  1/1500  km. 
Figure  6-2  shows  the  results  for  five  different 
distributions  of  data. 

Figure  6-2(a)  shows  results  for  both  observations 
collocated  with  the  gridpoint.  Both  observations  receive 
the  same  weight,  but  this  weight  is  about  20%  less  than 
the  weight  given  to  a  single  bathy  observation  made 
at  the  gridpoint  for  the  combination  of  parameters 
considered  (see  Table  6-2).  This  weight  reduction  is 
because  there  is  redundancy  in  the  two  observations 
as  a  result  of  their  collocation.  Nevertheless,  the  error 
achieved  for  this  case  is  the  minimum  of  all  the  test 
cases  as  the  observations  are  the  most  highly  correlated 
with  the  gridpoint. 


Figure  6-2(b)  shows  results  for  observation  1 
displaced  300  km  to  the  left  of  the  gridpoint  and 
observation  2  displaced  300  km  to  the  right  of  the  grid- 
point.  Because  of  the  symmetrical  distribution  of  the 
data,  both  observations  receive  the  same  weight. 
The  error  is  higher  than  that  of  Case  2(a)  because  the 
observations  are  less  correlated  with  (i.e.,  farther  away 
from)  the  gridpoint  than  before. 

Figure  6-2(c)  shows  results  for  observation  1 
displaced  600  km  to  the  left  of  the  gridpoint  and 
observation  2  displaced  300  km  to  the  right  of  the  grid- 
point.  Because  it  is  more  highly  correlated  with  (i.e., 
closer  to)  the  gridpoint  than  observation  2,  observation 
1  receives  a  higher  weight.  The  analysis  error  for 
Case  2(c)  is  higher  than  that  of  Case  2(b)  because 
observation  1  is  farther  away  from  the  gridpoint  in 
Case  2(c). 

Figure  6-2(d)  shows  results  for  both  observations 
displaced  300  km  to  the  left  of  the  gridpoint.  Both 
observations  receive  the  same  weight,  but  this  weight 
is  lower  and  the  analysis  error  is  higher  than  in 
Case  2(b).  These  results  reflect  the  higher  redundancy 
of  the  observations  caused  by  their  collocation. 

Finally,  Figure  6-2(e)  shows  results  for  observations 
1  and  2  displaced  600  and  300  km  to  the  left  of  the 
gridpoint,  respectively.  Note  that  the  weight  given  to 
observation  1  is  approximately  30%  less  than  in 


Case  2(c),  since  observation  2  falls  between 
observation  1  and  the  gridpoint.  Note  also  that  the 
analysis  error  is  higher  in  Case  2(e)  than  in  Case  2(c) 
because  both  observations  are  on  the  same  side  of  the 
gridpoint  (and  thus  more  redundant)  in  Case  2(e). 

In  general,  the  results  of  Figure  6-2  illustrate  the 
nonlinear  aspect  of  the  Oi  data  assimilation  technique. 
That  is,  the  weights  assigned  to  observations  depend 
strongly  on  the  distribution  of  the  observations  relative 
to  each  other,  as  well  as  relative  to  the  analysis 
gridpoint. 

7.0  Horizontal  Analysis  Grids 

OTIS  can  function  on  the  standard  FNOC  63  x  63 
northern  hemisphere  and  southern  hemisphere  polar 
stereographic  grids.  These  grids  cover  the  entire 
northern  and  southern  hemispheres  with  mesh  lengths 
varying  from  about  380  km  in  high  latitudes  to  about 
200  km  in  the  tropics.  In  addition,  it  can  run  on 
relocatable  63  x  63  or  125  x  125  regional  polar  stereo¬ 
graphic  grids  of  arbitrary  location,  extent,  and 
orientation.  Finally,  OTIS  can  function  on  a 
relocatable  rectangular  latitude/longitude  grid  covering 
an  arbitrary  domain  with  dimensions  up  to  80  x  180. 
Grid  spacings  for  regional  applications  of  OTIS  are 
typically  20-40  km. 
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Figure  6-2.  Analysis  error  ok  and  weights  given  to  two  bathy  SST  observations  a,  and  a2>  TOPS  ft,  and  climatology  y,  for  observa¬ 
tions  made  at  analysis  time  following  an  extended  period  of  no  data.  The  observational  locations  are  indicated  by  the  o's  and  the 
gridpoint  location  is  indicated  by  the  X.  The  scale  along  the  bottom  of  each  figure  shows  distance  in  kilometers.  All  results  assume 
A’’  =  0.5,  Gt  =  0.2°C,  a\  =  1°C,  and  Ak  =  1/1500  km. 
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8.0  Supporting  Climatologies 

OTIS  can  utilize  either  the  FNOC  hemispheric 
climatology  (Mendenhall  et  al.,  1980)  or  the 
Generalized  Digital  Environmental  Model  climatology 
(GDEM;  Davis  et  al.,  1986)  to  supply  Tji.  The  FNOC 
climatology  is  defined  by  month  on  the  standard 
FNOC  northern  hemisphere  and  southern  hemisphere 
63  x  63  polar  stereographic  grids.  The  GDEM  clima¬ 
tology  is  defined  by  season  on  a  0.5°  latitude/longitude 
grid  encompassing  the  north  Atlantic,  the  north  Pacific 
and  the  Mediterranean.  Both  climatologies  exist  in  disk- 
resident  form  on  the  FNOC  computer  system  and  are 
automatically  interpolated  in  time  and  space  to  the 
appropriate  analysis  time  and  grid  by  OTIS. 

9.0  Ocean  Feature  Models 

To  account  for  the  subsurface  thermal  structure  of 
ocean  fronts  and  eddies,  ocean  feature  models  can  be 
automatically  embedded  in  the  OTIS  first-guess  field 
at  known  front  and  eddy  locations  obtained  from 
satellite  imagery.  These  models  are  basically  canonical 
three-dimensional  ocean  front  and  eddy  structures 
based  on  historical  cruise  data  and  dynamical  insight 
(see  Spall,  1986;  Bennett  et  al.,  1989).  This  feature- 
model  embedding  process  produces  subsurface  thermal 
structures  in  the  vicinity  of  fronts  and  eddies  that  are 
dynamically  consistent  with  surface  currents  and 
substantially  more  accurate  than  climatology,  even  in 
the  absence  of  any  subsurface  data.  In  conjunction  with 
the  ocean  feature  models,  the  data  decorrelation  option 
in  OTIS  is  utilized  to  confine  the  influence  of 
observations  to  the  water  masses  in  which  they  are 
made.  Thus,  data  are  not  “smeared”  laterally  across 
fronts  and  the  boundaries  of  eddies,  resulting  in  a 
sharper  depiction  of  these  features. 

10.0  Example  Output 

In  this  section  we  present  output  from  the  global- 
scale  implementation  of  OTIS.  The  analysis  shown  here 
is  generated  on  the  FNOC  63  x  63  northern  and 
southern  hemispheric  polar  stereographic  grids.  The 
graphics  are  produced  by  interpolating  the  fields  to  a 
2.5°  spherical  grid  and  then  contouring.  Because  of 
the  coarse  resolution,  mesoscale  fronts  and  eddies  are 
not  represented  in  these  products. 

Figures  10-1  and  10-2  show  SST  and  SST  anomaly 
(analyzed  temperature  minus  climatology)  valid  on 
4  February  1988.  The  contour  interval  for  SST  is  1°C. 
The  contour  interval  for  SST  anomaly  is  0.5°C; 
negative  anomalies  are  shown  with  a  solid  line,  and 
positive  anomalies  are  shown  with  a  dashed  line.  As 
indicated  by  Figure  10-2,  prominent  large-scale 
anomalies  are  present  throughout  the  world’s  oceans. 
Of  note  are  (1)  the  two  cold  anomalies  in  the 


midlatitude  north  Pacific  larger  than  -2.5 °C  and 
centered  at  42°N,  172°E  and  36°N,  153°E;  (2)  the 
warm  anomaly  in  the  western  tropical  north  Pacific 
larger  than  2.0°C  and  centered  at  25°N,  152°E;  and 
(3)  the  warm  anomaly  in  the  midlatitude  south  Atlantic 
larger  than  3.0°  and  centered  at  32°S,  7°W. 

Figures  10-3  and  10-4  show  temperature  and 
temperature  anomaly  at  the  200-m  depth  from  the 
4  February  1988  analysis.  As  before,  the  contour 
interval  for  temperature  is  1  °C,  the  contour  interval  for 
temperature  anomaly  is  0.5°C,  negative  anomaly 
contours  are  displayed  with  a  solid  line,  and  positive 
anomaly  contours  are  displayed  with  a  dashed  line. 
From  Figure  10-4,  large-scale  anomalies  at  the  200-m 
depth  are  prominent  in  the  Pacific.  Of  particular  note 
are  (1)  the  cold  anomaly  in  the  western  tropical  Pacific 
larger  than  -3.5°C  and  centered  at  7°N,  155°E;  (2)  the 
two  warm  anomalies  in  the  eastern  tropical  Pacific 
larger  than  1.5°C  and  centered  at  20°N,  139°W  and 
17°N,  155°W;  and  (3)  the  warm  anomaly  in  the 
tropical  south  Pacific  larger  than  3.0°C  and  centered 
at  12°S,  142°W. 

11.0  Validation 

As  increasingly  sophisticated  ocean  thermal  models 
become  available  for  operational  use,  it  is  imperative  to 
have  the  means  of  rigorously  and  objectively  establish¬ 
ing  their  relative  skill.  This  validation  requires 
consistent  comparison  of  the  models  against  the  same 
independent,  accurate  ocean  thermal  data.  A  system  to 
carry  out  this  process  automatically  and  routinely  using 
bathy  data  is  in  place  at  FNOC  (SASC,  1984,  1985; 
Clancy  et  al.,  1986;  Clancy  et  al.,  in  press).  This  system 
is  used  for  validating  OTIS  and  other  ocean  thermal 
models  at  FNOC  and  provides  a  detailed  and  complete 
archive  of  model  performance.  Results  of  these  valida¬ 
tions  are  published  in  the  FNOC  Quarterly  Model 
Performance  Summaries,  and  are  used  to  ensure  that 
new  models  or  techniques  implemented  operationally 
are  indeed  more  skillful  than  the  old  ones  they  replace. 

The  thermal  fields  from  OTIS,  TOPS,  and  other 
models  are  interpolated  once  per  day  to  the  locations 
of  all  bathy  observations  received  at  FNOC  during 
the  previous  24  hours,  but  not  yet  assimilated  into  the 
models  (making  the  observations  independent  of 
the  models).  Each  bathy  observation  is  then  vertically 
interpolated  to  a  standard  grid  between  the  surface  and 
300  m,  and  spurious  observations  are  automatically 
discarded  by  a  bathy  error  detection  algorithm  (see 
SASC,  1985). 

A  full  suite  of  error  statistics  for  temperatures  at 
depths  of  0,  12.5,  25,  50,  75,  100,  150,  200,  and  300  m 
are  calculated  from  the  remaining  bathy  observations 
in  various  fixed  geographical  regions  (e.g.,  Eastern 
Pacific,  Western  Atlantic,  etc.).  These  statistics  include 
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Figure  10-  1 .  Global-scale  SST  from  OTIS  valid  on  4  February  1988.  The  contour  interval  is  I°C. 
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Figure  10-4.  Global-scale  temperature  anomaly  at  200  m  depth  (analyzed  temperature  minus  climatology)  from  OTIS  valid  on  4  February  1988.  The  contour  interval 
is  0.5  C  and  the  decimal  points  are  suppressed  in  the  labeling.  Thus,  for  example,  a  contour  label  of  25  implies  2.5  C.  Negative  contours  are  shown  with  a  solid  line 
and  positive  contours  arc  shown  with  a  dashed  line. 


the  root-mean-square  (RMS)  error,  the  mean  error, 
and  the  systematic  RMS  error. 

The  RMS  error  is  the  square  root  of  the  mean  of 
the  squared  differences  between  the  observed 
temperatures  at  various  depths  and  the  OTIS  analyzed 
temperatures.  Because  the  temperature  differences  are 
squared,  errors  of  opposite  sign  do  not  cancel  in  this 
statistic,  and  it  thus  supplies  information  on  the 
model’s  skill  at  predicting  the  pattern  of  the  thermal 
field.  Generally,  RMS  error  tends  to  be  contaminated 
by  observational  noise,  due  partly  to  instrumental 
errors  that  are  not  screened  out  by  the  error-checking 
algorithms  and  partly  to  subgrid-scale  errors  associated 
with  thermal  features  of  spatial  extent  too  small  to  be 
resolved  by  the  model  grid.  If  the  observational  noise 
is  large,  then  it  can  obscure  the  true  difference  in  the 
performance  of  models. 

The  mean  error  is  simply  the  average  of  the  predicted 
minus  observed  temperatures  at  various  depths.  It  gives 
an  indication  of  bias  in  the  predicted  thermal  fields 
over  the  entire  verification  domain,  but  tends  to 
average  out  model  errors  characterized  by  spatially 
varying  signs  (e.g.,  too  warm  in  part  of  the  domain 
but  too  cold  in  another).  Thus,  unlike  the  RMS  error, 
the  mean  error  supplies  no  information  on  how  well 
the  model  predicts  the  pattern  of  the  thermal  field. 
Observational  error  (which  has  a  zero  mean)  is 
suppressed  in  the  mean  error  statistic. 

The  systematic  RMS  error  is  defined  as  the  square 
root  of  the  mean  of  the  squared  differences  between 
the  observed  temperatures  and  the  least-squares 
regression  line  through  the  scatter  of  predicted  versus 
observed  temperatures  at  various  depths.  Like  the  RMS 
error,  this  statistic  provides  information  on  the  model’s 
ability  to  represent  the  pattern  of  the  thermal  field. 
Unlike  RMS  error,  however,  it  tends  to  suppress  the 
effects  of  observational  noise  in  the  verification  data; 
it  thus  allows  for  a  more  valid  intercomparison  of 
model  performance  than  the  RMS  error  (see  Willmott, 
1981,  1982). 

These  statistics  are  calculated  from  running  30-day 
and  90-day  histories  of  the  predicted  minus  observed 
temperatures,  and  are  output  in  graphical  form.  The 
resulting  plots  show  profiles  of  the  error  statistics  from 
the  surface  down  to  the  300-m  depth.  Figure  11-1  is 
an  example  that  shows  RMS,  mean,  and  systematic 
RMS  error  for  Northern  Hemisphere  (NHEM)  OTIS 
and  climatology  in  the  eastern  North  Pacific  (30-70°N, 
120-175°W)  for  the  month  of  January  1988  from  com¬ 
parison  against  406  independent  bathy  observations. 


OTIS  clearly  exhibits  substantial  improvements  over 
climatology  in  all  three  sets  of  statistics. 

12.0  Summary  and  Future  Developments 

OTIS  is  an  ocean  thermal  analysis  system  designed 
for  operational  use  at  FNOC.  It  is  based  on  the 
Optimum  Interpolation  data  assimilation  technique  and 
functions  in  an  analysis-forecast-analysis  data 
assimilation  cycle  with  TOPS.  OTIS  provides  a  rigorous 
framework  for  combining  real-time  data  (e.g.,  ship, 
buoy,  bathy,  and  satellite  observations),  predictions 
from  an  ocean  forecast  model  (e.g.,  TOPS,  OCEANS, 
etc.),  and  ocean  thermal  climatology  (e.g.,  GDEM)  to 
produce  a  statistically  optimum  representation  of  ocean 
thermal  structure.  OTIS  outputs  the  thermal  field  on 
a  vertically  stretched  grid  with  either  17  fixed  levels  in 
the  upper  400  m  (shallow  analysis  mode)  or  34  fixed 
levels  in  the  upper  5000  m  (deep  analysis  mode),  and 
can  function  on  hemispheric  or  relocatable  regional 
grids  that  employ  either  polar  stereographic  or  spherical 
projections.  Figure  12-1  illustrates  the  central  role  OTIS 
plays  in  linking  various  data  bases  and  models  to 
produce  a  comprehensive  ocean  thermal  analysis  and 
prediction  system. 

OTIS  Version  1  became  operational  at  FNOC  on  the 
NHEM  and  SHEM  63  x  63  polar  stereographic  grids 
in  July  1988  (see  Clancy  et  al.,  in  press).  It  was  upgraded 
to  function  on  the  FNOC  NHEM  and  SHEM  125  x  125 
polar  stereographic  grids  in  July  1989,  thereby  doubling 
the  spatial  resolution  of  FNOC’s  global-scale  represen¬ 
tation  of  ocean  thermal  structure.  OTIS  Version  2, 
designed  only  for  regional-scale  eddy-resolving  applica¬ 
tions  and  relying  heavily  on  use  of  ocean  feature 
models,  became  operational  at  FNOC  for  the  Gulf 
Stream  region  in  June  1989.  OTIS  Version  3,  which 
employs  a  more  sophisticated  use  of  ocean  feature 
models  and  supporting  data  bases,  is  currently  being 
implemented  to  run  on  the  CYBER  205  and,  eventually, 
the  Large  Scale  Computer  (LSC)  at  FNOC. 

A  future  development  will  involve  the  use  of  satellite 
altimetry  data  to  construct  synthetic  bathy  profiles  for 
use  by  OTIS.  This  data  may  prove  to  be  an  important 
source  of  thermal  information  below  the  400-m  depth. 
In  addition,  the  ocean  feature  models  and  supporting 
data  bases  used  by  OTIS  to  add  definition  to  the  sub¬ 
surface  thermal  field  associated  with  fronts  and  eddies 
will  be  improved  and  extended  to  all  regions  of  Navy 
interest.  Finally,  the  statistical  parameters  required  by 
OTIS  will  be  refined  as  experience  with  the  system 
grows. 
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Figure  11-1.  (a)  RMS,  ( b )  mean,  and  (c) 
systematic  RMS  temperature  errors  for 
NHEM  OTIS  (solid  tine)  and  climatology 
(dashed  line)  for  the  eastern  North  Pacific 
( 30-70° N,  120-175°W)  during  the  month  of 
January  1988.  The  statistics  are  based  on 
comparison  against  406  independent  bathy 
observations  and  the  error  bars  indicate 
approximate  90%  confidence  limits. 
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Appendix 

Basic  Theory  of  Optimum  Interpolation 


Here  we  derive  equation  (2-2)  of  Section  2,  which  is  the  fundamental  equation  of  OI  as 
implemented  in  OTIS.  Let  fl  represent  the  z'th  estimate  of  a  field  to  be  mapped  to  a  discrete  grid. 
Estimates  are  obtained  from  observations  and  the  predictions  from  a  numerical  model,  and  we 
can  write 


f,  =  F,  +  E,  (A-l) 

where  Ft  is  the  “true”  representative  value  of  the  field  on  the  scale  resolved  by  the  grid  and  Ej 
is  the  error  in  the  /th  estimate.  If  /,  is  an  observation,  then  £(  is  the  observational  error  due  to 
the  combined  effects  of  both  instrumental  and  aliasing  (i.e. ,  subgrid-scale)  errors.  If  ft  is  the 
predicted  field  from  a  numerical  model  functioning  on  the  same  grid  as  the  target  analysis,  then 
E:  is  simply  the  grid-scale  error  in  the  model  prediction. 

The  deviation  of  the  mapped  field  from  a  first-guess  field  is  represented  as  a  weighted  combination 
of  the  estimated  departures  from  the  first  guess.  This  representation  is  subject  to  error,  and  we 
can  thus  write 

f;  =  si'f;p,  +  /*  (a-2) 

i  =  i 

where  Fj  is  the  true  representative  value  at  gridpoint  k,  pk  is  the  weight  given  to  the  z'th  estimate, 
N  +  1  is  the  number  of  estimates  assimilated  at  the  gridpoint,  and  Ik  is  the  analysis  error  at  the 
gridpoint.  The  primes  denote  departure  from  the  first-guess  field,  which  is  assumed  to  be  the  long¬ 
term  climatology  and  thus  an  unbiased  estimator. 

From  equations  (A-l)  and  (A-2),  the  mean-square  analysis  error  at  gridpoint  k  over  an  ensemble 
of  analysis  realizations  is  given  by 


r'iV  +  f,)p,-f;)2 


i  - 1 


(A-3) 


which  can  be  expanded  to 


K  = 


X '  1 '  p,Pj  (f;  f;  +  6ue;  e;  >  -  2  "i  'P,  f;  f’  + 

1  =  1  ;  =  1  i=l 


f;f; 


(A-4) 


where  the  overbar  denotes  an  ensemble  mean  over  many  analysis  realizations,  6tJ  is  the  delta 
function  (dJy  =  1  for  z  =  y;  dly  =  0  for  z  ¥>  j),  and  we  have  made  the  usual  assumptions  that 
the  correlations  between  the  errors  of  the  estimates  and  the  true  deviations  from  the  first-guess 
field  of  climatology  (i.e.,  £/  Fj  ),  and  the  correlations  between  the  errors  of  the  estimates  (i.e., 
£/  Ej  for  z  #  j)  are  zero  (see  Bergman,  1978).  Dividing  (A-4)  by  the  variance  of  F '  (i.e„  Fj  Fj  ) 
yields 


=  I  I  PiPjblij  +  <5,A>-2  I  PM*.  +  1 


/v 


(A-5) 


1=1  7=1 


i  -  1 


where  V  is  the  normalized  mean-square  analysis  error, 


n„  = 


F; 


and 


£/£/ 


f;  £/ 


(A-6) 


(A-7) 


The  quantity  n  is  the  space-time  autocorrelation  function  between  locations  i  and  j  in  sampling 
space  for  the  resolvable  signal  £',  and  A,  is  the  noise-to-signal  ratio. 

To  obtain  the  optimum  values  of  the  weights  Fj  applied  to  the  estimates  of  the  target  field/' , 
the  normalized  mean-square  analysis  error  H*  is  minimized  by  setting 


(A-8) 


From  (A-5),  this  yields 

I  (*1„  +  =  1*,  O'  =  1,  2,  3  ...  N+  1). 

j  -■  i 

Now,  in  the  context  of  OTIS,  the  first  N  estimates  are  observations  (i.e.,  ship,  bathy,  and  satellite 
data)  and  the  N  +  1th  estimate  is  a  prediction  from  the  TOPS  model.  Using  the  subscript  k  to 
denote  the  Aih  analysis  gridpoint  and  defining  the  weights  applied  to  the  observations  as  akJ, 
the  weight  applied  to  TOPS  as  flk,  the  noise-to-signal  ratio  for  the  observations  as  A(°,  and  the 
noise-to-signal  ratio  for  TOPS  as  Xk,  equation  (A-9)  reduces  to 


(A-9) 


i  ( n,j  +  <5,jK)akj  +  ( n ik  +  6,kXpk)pk  =  rj*, 

i  --  i 

for  i  =  1,  2  ...  ,  N,  (N  +  1)  =  k 
which  is  precisely  equation  (2-2)  of  Section  2. 


(A- 10) 
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